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We look at two possible routes to classical behavior for the discrete quantum random walk on the 
line: decoherence in the quantum "coin" which drives the walk, or the use of higher-dimensional 
coins to dilute the effects of interference. We use the position variance as an indicator of classical 
behavior, and find analytical expressions for this in the long-time limit; we see that the multicoin 
walk retains the "quantum" quadratic growth of the variance except in the limit of a new coin for 
every step, while the walk with decoherence exhibits "classical" linear growth of the variance even 
for weak decoherence. 
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INTRODUCTION 

Considerable work has been done recently on quantum 
random walks, which arc unitary (and hence reversible) 
systems designed as analogues to the usual classical case. 
Two approaches have been taken to the problem: con- 
tinuous 0] and discrete |^, ^, |^ walks. This paper is 
exclusively concerned with the discrete walk on the in- 
finite line. In this case, we introduce an extra "coin" 
degree of freedom into the system. As in the classical 
random walk, the outcome of a "coin flip" determines 
the direction that the particle moves. In the quantum 
case, however, both the flip of the coin and the condi- 
tional motion of the particle are unitary transformations. 
Different possible classical paths can therefore interfere. 

For the classical walk, p{x, t) has the form of a bino- 
mial distribution, with a variance — which grows 
linearly with time. The variance in the quantum walk, by 
contrast, grows quadratically with time [||; and the dis- 
tribution p{x,t) has a complicated, oscillatory form [Q. 
Both of these are effects of interference between the pos- 
sible paths of the particle. There are two obvious ways to 
regain the classical behavior. If the quantum coin is mea- 
sured at every step, then the record of the measurement 
outcomes singles out a particular classical path. By aver- 
aging over all possible measurement records, one recovers 
the usual classical behavior. 

Alternatively, rather than measuring the coin every 
time, one could replace it with a new quantum coin for 
each flip. After a time t one would have accumulated 
t coins, all of them entangled with the position of the 
particle. By measuring them, one could reconstruct an 
unique classical path; averaging over the outcomes would 
once again produce the classical result. 

These two approaches, which are equivalent in the clas- 
sical limit, give two different routes from quantum to 
classical. We might increase the number of coins used 
to generate the walk, cycling among M different coins, 
in the limit using a new coin at each step. Or we might 



weakly measure the coin after each step, reaching the clas- 
sical limit with strong, projective measurements. This is 
equivalent to having a coin which is subject to decoher- 
ence. 

In this paper we contrast these two approaches, using 
the functional dependence of the variance on time as an 
indicator of "classical" vs. "quantum" behavior. In the 
presence even of very weak decoherence, the variance of 
the quantum walk grows linearly with t at long times; 
while even using a large number of coins, the the variance 
of the unitary walk grows quadratically. 

Let us consider a fairly general quantum walk on the 
line. The particle degree of freedom has a basis of posi- 
tion eigenstates {js^)}, x\x) = x\x), where x can be any 
integer. We will assume that the particle begins the walk 
at the origin, in state |0). The walk is driven by a sep- 
arate "coin" degree of freedom: a D-dimensional system 
with an initial state |$o)- Let 'Pr,'Pl be two orthogonal 
projectors on the Hilbert space of the "coin," such that 
Vk + 'Pl = I- These represent the two possible outcomes 
of the coin flip. Right or Left. We assume that the coin 
is unbiased, meaning TtVr = TiVi, = D/2. We also de- 
fine a unitary transformation U which "flips" the coin. 
One step of the quantum walk is given by the unitary 
operator 

E=(^S®Vk + S^® Vl) (/ C>) , (1) 

where S, are shift operators on the particle position, 
S\x) = \x + l), S'flx) = \x-l). The full initial state of 
the system (particle and coin) is I'^'o) = |0) I'l'o)- 
We can identify the eigenvectors |fc) of S, S\ 

k> = r ^e-^'^\k) , (2) 

with eigenvalues 

5|A:) =e-*'=|fc) , S^\k) = e+'^'lk) . (3) 
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In the k basis, the evolution operator becomes 

= |fc)®;7fe|$) , 



(4) 



where Uk is also a unitary operator. 

We now generalize to allow for decoherence. Suppose 
that before each unitary "flip" of the coin, a completely 
positive and unital map is performed on the coin. This 
map is given by a set of operators {v4„} on the coin degree 
of freedom which satisfy 



(5) 



A density operator x for the coin degree of freedom is 
transformed 



(6) 



We combine this with the unitary evolution to define an 
evolution superoperator 



(7) 



Note that for the diagonal case k = k', this superoperator 
is also unital, and hence preserves the identity. 
The initial state is 

/Hk r dk' 
_y _|fc)(A:'|®|$o)(<i>o|. (8) 

Let the quantum random walk proceed for t steps. Then 
the state evolves to 

/§/ ^l^>(^'l®4.'l<i>o>($o|. (9) 
' to reach a point x at time t 

x\k){k'\x)i:T{cl^,\^o){M] 



pt 



From this, the probability to reach a point x at time t is 

dk r dk' 
2n J 27r ' 

dk r dk' 



p{x,t) 



2tt / 2tt 



""''-''^Tr{4,,|$o)($o|KlO) 



Eq. ( p^ ) for p{x,t) will be difficult to evaluate analyt- 
ically; hence, we restrict our interest to the moments of 
this distribution. 

X 

dk f dk' 



X 

xTr{4,,|<i>o)(<i>o|} 



-ix{k—k') 



(11) 



We invert the order of operations and do the x sum first, 
yielding 



xTr{4,,|ci>o)(<f>o|} 



dk / dk'6'^"'\k-k') 



(12) 



where 5^'^\k — k') is the mth derivative of the delta func- 
tion. We can then integrate this by parts. 
In integrating ( p^ ) we will need 

^TrjA-fc^o} = -zTr{zAfc,o} , 
= -zTr{(£fcfc,0)z} , 



dk 

Z = Vn-PL 



(13) 



Making use of ([l3|), when we carry out the integration 
by parts for the first moment we get 



^^Tr{z4|$o)(<i>o|} 



dk 
2 



(14) 



where we've introduced the simplified notation Ck = Ckk- 
We can carry out a similar integration by parts to get 
the second moment: 



/i[EETr{^/:r'(^4>o)(<i>oi)} 

+ E E Tr{z/:p'((£i.>o)($o|)^)} 



Let us for the moment specialize to the unitary case, 
so CkP = UkpUl- In this case, we can expand |<I>o) in 
terms of the eigenvectors of Uk- 



kll 



(16) 



Assume, for the moment, that Uk is nondegenerate, so 
the 9ki are all distinct. If we plug these expressions into 
( |l^ ) and ([TbI), we notice that most of the terms will be 
oscillatory; at long times i, they will average to zero. The 
only terms that survive will be diagonal in I: 



J2 h^('^o\iUkyziuly\<Po) 



-t 



D 



^Y.\^ki?{Mz\M 



1=1 

-|-oscillatory terms . 
Similarly, for the second moment 

D 



(17) 



f 



-0{t) + oscillatory terms . 



(18) 



So in the long-time limit, the variance will always grow 
quadratically in time for a unitary coin of finite dimen- 
sion. If Uk is degenerate, the expressions ( pj| ) and ( p^ ) 
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FIG. 1: The variance (a;^) — {x)'^ for the unitary walk with 
multiple coins, Af=l-5. The coins all begin in state |_R). The 
variance grows quadratically in time, in contrast to the linear 
growth in the classical random walk. 



will have to be modified to include appropriate cross- 
terms; but this will not change the functional dependence 
on t. 

The usual case considered in the literature has taken 
the coin to be a simple two-level system, and the "flip" 
operator U to be the usual Hadamard transformation H: 
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H\R) = —{\R) + \L)) , 
H\L) = 1={\R)-\L)) . 



(19) 



The projectors are Vn = \R){R\, = \L){L\. The walk 
on the line in this case has been exactly solved by Nayak 
and Vishwanath and agrees with the expression dl^ ) 
given above. We have also considered the case of a walk 
driven by M coins, flipped cyclically In this case, 
and ^ can be solved analytically; at long times, 
the variance goes like 

(£2) _ {if ^ f |^^^.^^±1/^ j + 0{t) + osc. terms . 

(20) 

We compare this to the results of numerical simulations 
in figure 1, finding the agreement to be excellent. Note 
that the quadratic growth of the variance doesn't vanish 
even in the limit of large M . Only with a new coin for 
every step {M = t) do we recover the classical behavior. 
(Note that this doesn't contradict the result ([2C|), since 
that is only strictly valid for t ^ M.) 

Let's now allow for decoherence. Because the superop- 
erator Ck is unital it preserves the identity Ckl — I', its 
largest eigenvalue is 1. We explicitly assume that all the 
other eigenvalues of Ck obey |A| < 1. It is handy then to 



separate out the traceless part of the initial state 
po = + (|$o)($o| - I/D) = i/D + xo . 
Plugging this into the equation (|lj) we get 



(21) 



3 = i 

dk 



- J §Tr{z(l-A)-i(£,-4+i)xo}(22) 

Note that (1 — £fe)~^ acts only on the traceless operator 
XOj so it is well-defined. If the coin is unbiased (as we 
assumed), then TrZ = and the first term vanishes. At 
long times the JO-I^^ term will decay away. So in the long 
time limit, the first moment of the walk with a decoherent 
coin will tend to a constant. 

The second moment is more complicated in detail, but 
similar in spirit. Again separating out the traceless part 
Xo, we get 



dk y--^ 

9^ 2^ 



Tr 



+Tr{zc{-^{Z{Cixo) + {J^ixo)Z)} 



dk X ^ 

2^^ 



1 + {2/D)TTiZ{l- C 



x(A-4--'+i)Z 



+Tr<{ Z{l-£k)-\Ck-£l-'+') 



x{ZiCixo) + iCixo)Z) 



(23) 



It is not difficult to see that the last term will tend to- 
wards a constant for large t, while the first two will grow 
linearly. So we get the approximate expression at long 
times 



1 + ^y" dkTT{z{l^Ck)-'CkZ}^ 



-const. , 



(24) 



which is linear in t. 

If we specialize to the case of a single two-level coin 
undergoing the Hadamard evolution, we can find exact 
solutions to (p5). We first pick a particular form for the 
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FIG. 2: (x^) — (i)^ vs. t for the quantum random walk with 
decoherence, for d = 7r/16, tt/S, 37r/16, 7r/4. For all cases the 
coin began in the initial state |_R). Note that the variance 
goes asymptotically to a linear growth at long times which 
matches our analytical estimate; the rate of growth goes to 
one with increasing decoherence, matching the classical case 
at 6) = 7r/4. 



decoherence. The most convenient numericahy is 'pure 
dephasing 



1 



(25) 



For 9 — this reverts to the unitary case, while for 6 — 
7r/4 there is complete decoherence at each step. At long 
times we have IItI 



~ {x)lK.t {co\? 29 + csc^ 26 



const. 



(26) 



which goes to 1 as ^ 7r/4, and diverges as 6* — + 0, when 
the long-time approximation breaks down. In figure 2, 
we compare this result to the output of numerical simu- 
lations, once more finding excellent agreement. Note that 
while the variance grows linearly, as in the classical case, 
it grows faster than the classical case. The reflects the 
persistent effect of interference, which causes the particle 
to continue to drift in a particular direction. 

From the time-dependence of the moments, one might 
reasonably claim that the multicoin system remains 
"quantum" even in the limit of very large numbers of 
coins, while the decoherent system remains "classical" 
even in the limit of very weak noise. Moreover, from 



equations and ( |24| ) this behavior seems to be generic, 
independent of the particular model chosen. 

We should emphasize that in this paper we have only 
altered the coin degree of freedom. One might naturally 
consider modifications of the evolution of the particle as 
well, such as allowing decoherence of the position as well 
as the coin. There have been numerical studies of this, 
and of the transition to classical behavior that results . 
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